draft version 



Evolution of cross-correlation and time lag of 4U 1735-44 along 

the branches 

Ya-Juan Lei 1 , Hao-Tong Zhang 1 , Cheng-Min Zhang 1 , Jin-Lu Qu 2 , Hai-Long Yuan 1 , Yi-Qiao 
Dong 1 , Yong-Heng Zhao 1 ,De-Hua Wang 3 , Hong-Xing Yin 4 , Li-Ming Song 2 

ABSTRACT 

We analyze the cross-correlation function between the soft and hard X-rays 
of atoll source 4U 1735-44 with RXTE data, and find the anti-correlated soft 
and hard time lags of about hecto-second. On the island state, the observations 
do not show any obvious correlations, and most observations of banana branch 
show positive correlation. However, the anti-correlations are detected on the 
upper banana branch. These results are different from those of Z sources (Cyg 
X-2, GX 5-1), where the anti-correlation is detected in the low luminosity states, 
then the lag timescales of both this atoll and Z sources are found to be similar, 
at the magnitude of several tens to hundreds of seconds. As a comparison, it is 
noted that the anti-correlated lags of thousand-second have been reported from 
the several black hole candidates in their intermediate states. Finally, we compare 
the correspondent results of atoll source 4U 1735-44 with those observed in black 
hole candidates and Z sources, and the possible origins of the anti-correlated time 
lags are discussed. 
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Introduction 



The timing and spectral analysis of the X-ray emission of the compact objects, neutron 
stars (NSs) or black holes (BHs 
ray binaries (LMXBs) (see, e.g. 



— „ „„„„„ 

van der Klis 


2006; 


Liu et al. 


2007a) 



spectral properties, which are shown in the defined color-color diagram (CCD) or hardness- 
intensity diagram (HID), the NS LMXBs can be classifi ed into two subtypes, Z source s with 
high luminosity and atoll sources with low luminosity (IHasinger fc van der Klislll989l ). The 
main difference between atoll and Z sources might be the accretion rate and magnetic field, 
where Z source shares the higher accretion rate and slightly stronger magnetic field than 
those of atoll source (jZhangl 120071 ) . The atoll sources cover a wider range of luminosity. 



wher eas Z sources usually appear as the Eddington limited luminosities (e.g., Ford et al. 
20001 ) . Z (atoll) source traces out a Z-shape track (island + banana pattern) within a few 
days (days to weeks). In CCD, the Z-track consists of three main branches, horizontal branch 



(HB) , normal branch (NB) and flaring branch (FB) ( IHasinger fc van der Klislll989t IHasinger 
1990). The source evolves continuously along the Z-track, with each position related to the 
different mass accretion rate. The accretion rate increases from HB to NB and reaches a 



maximum at the right end of the FB (IHasinger et al.lll990l ) . However, recently, the results of 
studying the source XTE J1701-462 suggest that the motion along the Z branches might be 
corresponding to the roughly constant accretion rate ( iLin et al.ll2009l ; iHoman et al.ll2010l ). 



As for the CCD pattern of atoll source, it follows the similar correspondences to those 
of Z source, from isla nd state (IS) to th e lower left banana (LLB ), lower banana (LB) and 
upper banana (UB) (Ivan der Klisl l2000t lAltamirano et ajj|2008j ). In the atoll source, the 
count rates increase with the inferred accretion rate M, which increases from the IS to the 



2000. 


2006; 


Belloni et al. 


2007) 



account of the low sensitivity at the low count rates. 

The spectral characteristics of Black hole X-ray binaries (BHXBs) usually exhibit five 
distinct states, quiescent state, low/hard state(LHS), inte r mediate state (IMS), high/soft 



state (HSS) and very high state (VHS) (IRemillardl l2005t iRemillard fc McClintockl 12006 
Bellonill2010[ ). The I MS is complex, who s e luminosity cou ld not be determined solely by the 
mass accretion rate ( Yu Sz Dolencel 2007 ; Yu Sz Yan 2009 ). Although the luminosity of VHS 
is much higher than that of IMS, both states share the similar spectral and timing behaviors, 
which can be take n as the same state that represents the transitions between LHS and HSS 
jDoneet al.ll2007j) . 



The timing properties of X-ray binary are related to the spectral states i.e., the positions 
of CCD. The cross-correlation function of the soft and hard X-ray photons can be used for 
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analyzing the relation of different energy bands, which is helpful for understanding the 
radiative and geometrical structure of the accretion disk for different spectral states. The 
cross-correlated soft and hard X-ray lags have been investigated for Z sources and BHXBs. 
For some BHXBs, the anti-correlation of a few hundred seconds a re often seen in the steep 



powe r-law (SPL) or intermediate states (jChoudhury et al.l 120051 ; ISriram et al.l 120071 12009 



20101 ) , and these timescales of the lags could be ascribed to the viscou s timescales of matter 



flows in the optically thick accretion disks (jChoudhury fc Raol 12004 ). Such anti-correlated 
lags (~ sev 



(Lei et al. 


2008; 


Sriram et al. 


2012) 



indicates that the phenomena of observed anti-correlations should not arise from the spurious 
nature, but the reality of physical processes. 

In some aspects, the atoll source and BH system share the similar spectral characteris- 
tics ( Idi Salvo et al.l 120061 ) . thus we select an atoll source, 4U 1735-44, to study its evolution 
of cross-correlation and time lag between the soft a nd hard X-rays. This source is classified 



as an atoll source with a relatively high luminosity (ISmale et al. 



distance of which is measured to be about 8.5 kpc ( iGalloway et al.l 120081 ). T he binary pa 



1986 



Seqn et al.lll997l) . the 



rame ters of 4U 1735-44 have been observed with th e orbital period of 4. 654 hr (jCorbet et al. 



19861 ) and its co mpanion mass o f 0.5 3±0.44 (Casares et al.l 120061 ). Its X-ray spectra 



are analyzed by iNg et al.l (120101 ) and iMueck et al.l (120111 ). with an iron emission line de- 
tected. In addition, the kHz QPOs are disc overed at the relatively lower mass accretion 



rates (IFord et al.l Il998l ; IWijnands et al.lll9981 ). In detail, we investigate whether the anti- 
correlation between the soft (2-3.3 keV) and hard (12-30 keV) X-rays exists in 4U 1735-44. 
As a further step, we compare the results of atoll with those of Z sources and BHXBs. The 
organization of the paper is described below. Section 2 is about the observations and data 
analysis of RXTE, and the results on the cross-correlation between the soft and hard X-rays 
are analyzed in section 3. The discussions and summary are written in section 4. 



2. Observation and Data Reduction 

We analyze all pointed observations of the Proportional Counter Array (PCA) on board 
the RXTE satellite. The PCA consists of 5 non-imaging, coaligned Xe multiwire proportional 
counter units (PCUs). In the work, only PCU2 data are adopted, which is the best calibrated 
unit and has the longest observational duration. After excluding the data segments less than 
2000 s, we extract the light curves from the data of standard 2 mode with bin size of 16 
s. With the XRONOS tool ll crosscor'\ we estimate the coefficient of the cross-correlation 
between the soft (2-3.3 keV) and the hard (12-30 keV) X-rays. 
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Similar to the previous publication about the cross-correlation for Z source Cyg X-2, 
we divide the cross -correlation res ults into three groups, the positive, ambiguous and anti- 
correlated (also see lLei et al.ll2008l ). We can obtain the time lag by fitting the anti-correlated 
part with an inverted Gaussian function. After subtracting the background, the light curves 
and anti-correlations between soft and hard X-rays, accompanying with the pivoting spectra 
from the data with different hardness ratios, are shown in Figure 1, where we define the hard 
and soft regions according to the hardness ratios. 

In order to study the source evolution on the CCD of 4U 1735-44, we define the soft 
and the hard colors as the count-r ate ratios 3.5-6.0 keV/2.0-3.5 keV and 9.7-16 keV/6.0-9.7 
keV (also see lO'Brien et al.ll2004l ). The light curves are extracted from the 2.0-3.5 keV, 
3.5-6.0 keV, 6.0-9.5 keV and 9.5-16 keV, respectively, with the background subtracted. By 
the correspo nding Crab values, we cor rect the affections of the different gain epochs and the 
gain change (Ivan Straaten et al.ll2005l ). Figure 2 shows the revised CCDs with the bin size of 
512 s. For spectral analysis, the PCA background subtraction is carried out using the latest 
versions of the appropriate background models, and a 0.6% systematic error is added to the 
spectra to account for the calibration uncertainties. As usual, the spectral fitting software 
XSPEC is used. 



3. Results 



3.1. The results of cross-correlations and their distribution on CCD 



With the data of RXTE, for the atoll source 4U 1735-44, we analyze the cross-correlation 
between the soft (2-3.3 keV) and hard (12-30 keV) X-rays. For each segment of the light 
curves, we analyze the cross-correlation. The ant i- correlations have been detected in the 
light curves of 12 observations. We apply the letter 'a' to mark the segment of the anti- 
correlation. The ranges of hard and soft time lags for the anti-correlation are from a few 
tens to hundred seconds. There are four regions in Figure 2 that shows the CCD of all 
analyzed observational data. The data in region T, 'II', 'III' and 'IV are mainly related to 
the IS, LLB, LB and UB, respectively. 

In Table 1, we list the observations with anti-correlation detected, the cross-correlation 
coefficients and delay times. From Table 1 and Figure 2, we notice that most of the anti- 
correlated observations are located at the UB, where show the 7 observations with anti- 
correlated hard X-ray lags, 3 observations with anti-correlated soft X-ray lags and 2 obser- 
vations with unobvious time lag. The delay times are simila r to those found in the Z sources 
Cyg X-2 and GX 5-1 (ILei et al.l l2008t ISriram et al.l 120121 ). There is no relation detected 
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between the delay times and the location in region 'IV of Figure 2. 



From Figure 2 and Table 2, we can see that, in the region T, ambiguous correlated 
data are dominated, then the great majority of positively correlated data dominate in the 
regions 'II' and 'III', and in the region 'IV, all the three correlations are found in the equal 
weight. Averagely, the correlation coefficients of the region 'II' are mostly less than those of 
region 'III'. Obviously, most anti-correlated data occur in regions 'IV. In addition, for the 
observations with the anti-correlations, the spectral evolution of the soft and hard regions 
of the light curves is shown in Figure 1. However, the energy ranges of the pivoting are not 
obviously evolutional along the locations of the CCD. There is not yet found the pivoting 
in the spectral evolution of p ositive-correlatio n be tween the sof t and hard X-rays, which is 
consistent with the results of iLei et al.l ( 120081 ) and ISriram et al.l ( 120121 ). 



3.2. Spectral variation 



It is known from the previous work tha t the position in t he CCD generally is well 
correlated to its spectral/timing state (see, e.g.. Ivan der Klisll2006[ ) . For studying the spectral 
evolution during the positive and anti-correlation between soft and hard X-rays, we analyze 
the spectra of the segments a, b, c and d (account rates decrease from a to d) of ObsID 91025- 
01-04-03, separately. For segment a, the cross- correlation between soft and hard X-rays shows 
anti-correlation, and for segment b the cross-correlation shows strong positive correlation, 
and, for segment c and d, the cross-correlations show obviously positive correlation but with 
a lower correlation coefficients than that of segment b. 

These spectra can be fitted b y the multi-color disk model (diskbb in XSPEC) plus 



black body (bbody in XSPEC) (e.g.. lWeng fc Zhang! 1201 ll ). A gaussian line with the center 



energy fixed at 6.6 keV and the width fixed at 0.5 keV is added to describe the iron line. 
The photo-electric absorption is modelled with wabs in XSPEC. Due to the absence of low 
energy observation of PCA, the absorption column is fixed at the value of Nn = 3.4 x 10 21 



cm 2 (IChristian &: Swank! Il997l ) 



Figure 3 and Table 3 show the fitting results, from segment a to d, as the X-ray flux 
decreasing, both the disk temperature and black body temperature decrease independently. 
We can see that, with the increasing of the flux (from d to a), the inner disk radius(oc 
the normalization) decreases, as the region of black body emission (a the normalization) 
increases. In addition, for segments a and d, both the soft and hard components dominate 
with the similar weights, the spectral p roperty of which is corresponding to the SPL states of 
BHXBs (also see lWeng &: Zhangj|201ll ). For segments b and c, the soft components dominate 
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the spectral character istic, which is similar to the HSS of BHXBs (also see lSriram et al.ll2010 
Weng fc Zhanifollh . 



4. Summary and Discussion 

We have detected the anti-correlation between soft (2-3.3 keV) and hard (12-30 keV) 
X-rays from the atoll source 4U 1735-44 in 12 RXTE observations. The results show that, 
the observations with anti-correlation are located at the UB, and the observations with 
positive correlation mostly appear in the LLB and LB. Most observations in the IS show the 
ambiguous correlation. The pivoting of the wide-band X-ray spectrum is also detected in all 
the observations with the anti-correlation detected. The energies of pivoting are generally 
lower than those detected in Z sources, and are not related to their location of CCD (see 
Table 1). 

The results from the atoll source 4U 1735-44 are different from those of Z sources 



(Lei et al. 


2008; 


Sriram et al. 


2012) 



locate at HB and upper NB, which correspond to "the low-hard state", however the anti- 
correlated observations detected in the atoll source 4U 1735-44 locate at UB which corre- 
sponds to high luminosity state. The anti-correlated time lag scales of Z and atoll sources 
are similar, and the observations with anti-correlated hard X-ray lags are more than those 
of soft X-ray lags for both Z and atoll sources. For Z source, if the accretion rate is the 
primary parameter to determine the location on CCD, the anti-correlation could be corre- 
sponding to the lower mass accretion. However, the m ass accretion rate may be constant 



for XTE J 170 1-462 during the Z-track (ILin et al.ll2009l ). Similar to Z sources, there is no 



obvious relation between the observed delays and corresponding positions on CCD for the 
atoll source 4U 1735-44. 

The anti-correlations have been detected in a few BHXBs with a few hundred to thou- 
sand seconds lags with the cross-correlation function of soft and hard X-rays, which suggests 



that the accretion disk could be truncated (jChoudhury et al.ll2005t ISriram et al.ll2007l . 12009 



2010l ). For two Z sources, Cyg X-2 and GX 5-1, the anti-correlation between soft and h ard X- 
rays favors the truncated accretion disk geometry (ILei et al.l 120081 ; ISriram et al.l 120121 ) . The 
truncated accretion disk model is ofte n used to explain the spectral and temporal features of 
the hard states observed in BHXBs JPone et_aLjl2007h however some observational results 



show that the disk may not be truncated ([Miller etal.1 12006c iRvkoff et al.l 120071). T h e iron 
line is corresponding to the innermost disk regions rtstellalll99oh . bone fc Gierlihskil (bood ) 
sugg est that truncated d i sk mo dels are consistent with the detected smeared iron emission 
line. iDone Diaz Trigol ( 120101 ) show that the iron line is instead narrow and consistent with 
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the truncated disk geometry in GX 339-4. For the atoll source 4U 1705-44, the similar profile 
of iron emission line is found in the high/soft states, which sugg ests that the accretio n disk 
could be also truncated in the high/soft states of atoll source (jdi Salvo et al.l 120091 ) . The 
anti-correlation detected in the high luminosity state should also imply that the truncated 
disk occurs there. 



For BHXBs, lEsin et al.l ( 119971 ) propose that an accretion flow around a BH consists of 
two zones with a transition radius, an inner advection-dominated accretion flow (ADAF) 
and an outer standard thin disk. In generally, in the LH state the disk is assumed to be 
truncated at a large radial distance, and in the HS state it becomes non-truncated. Another 
suggestion is that the accretion disk is truncated very close to the black hole in the IM 
state/SPL st ate from the spect ral and temporal results of the study of various black hole 



sources (e.g, 



Sriram et al.ll2010l ). The anti-correlation detected in the BHXBs implies that 
the ADAF condenses and expanses to the inner disk, and therefore the hard X-ray flux 
increases and soft X-ray flux decreases. The IM state is probably the spe ctral state where 



Mever et al. 


2007; 


Sriram et al. 


2010) 



The NS has a solid surface and magnetic field, which is different from those of BH. 
However, if the luminosity of NS LMXB increases above a critical luminosity, its magneto- 
sphere may go inside IS CO, w hich may account for the same disk evolutionary pattern for 
BHXBs IWeng &: Zhangl (120111 ). Our results show that the anti-correlation is detected in the 
UB which is corresponding to the max flux for the atoll source 4U 1735-44. According to the 
fitting results of spectra and X-ray luminosity, the LB could be corresponding to the HSS of 
black hole, and the UB is in accordance with the VHS. Although the different luminosity for 
IMS and VHS of BHXBs, both states share the similar spectral and timing behavior, and 
they are often tho ught as the same state and stand for the transitions between LHS and HSS 
(IDone et al.l 120071 ). The anti-correlations occur at SPL of BHXBs and at UB of NS system 
(atoll), so we can take the UB of atoll to be similar to SPL of BH system. 

Furthermore, we notice that the soft and hard X-rays are not correlated in the IS state 
for atoll source, which could be due to the fact that the accretion disk boundary is far 
away from the hot corona where emits the high energy photons. With the accretion rate 
increasing, the accretion disk moves inward with more soft seed photons produced, at the 
same time, the ADAF becomes be condensed and produces more hard photons, therefore 
the positive-correlation occurs, implying the increases of both the soft and hard X-rays. 
When the luminosity increases above a critical luminosity, the ADAF extends outside to the 
form of the corona, which expands quickly, so the soft X-rays decrease with the hard X-rays 
increasing. In the UB, the accretion disk extends down to ISCO and could be truncated 
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there, and the anti-correlations can be detected. Our spectral analysis suggests that with 
the increasing of the accretion rate, the inner radius of the accretion disk decreases, which is 
consist ent with BHXBs th at the accretion disk could be truncated very close to the compact 
object (jSriram et al.ll2010l ). 



The anti-correlated lags are detected and range from a few tens to h undred seconds 
in atoll source 4U 1735-44, which are similar to those found in Z sourc es (ILei et al.l 12 008; 
Sriram et al1l2012h but time lags are lower than those in BHXBs (e.g., IChoudhury fc Rad 



2004 



Sriram et al.l 120091 ) . Generally, the hard X-ray emission could come from the Comp- 
tonizing of the soft seed photons. For BHXBs, the soft seed photons could only come from 
the accretion disk. For NS LMXB systems, the soft seed photons could be produced from 
the surface of NS and/or the accretion disk. The Comptonizing might occur at a hot corona, 
a h ot flared inner disk, or even the boundary layer between the NS and accretion disk (also 



sec 



Popham fc Sunyaevl l200ll ; iQu et al.l 120011 ). The time scale of only Comptonization pro- 



cess is expected to be about <-, 1 second, therefore the observed several hundred-second lags 



could be the other processes (lHasingerlll987l ; iNowak et al.lll999uBottcher fc Lianglll999l ). We 



infer that this anti-correlated lag time scale of about several hundred seconds would be re- 
lated to the change of the disk structure. The different anti-correlated lag timescale implies 
that the disk inner radius occurs at various locations, which can also affect corona region. 



Sriram et al.l ( 120121 ) suggest that the observed lags could be due to the change in the size of 



the Comptonizing region. Some authors suggest that the lag time of several hundred sec- 
onds is comparable to 
processes change (e.g. 



Choudhurv & Rao 


2004; 


Lei et al. 


2008; 


Sriram et al. 


2012) 



ever, above model cannot explain the anti-correlated soft X-ray lag, which can be explained 
that a fluctuation produced from the innermost accretion disk could propagate to the outer 
accretion disk and affect the soft X-ray emission regions. The fluctuation of accretion disk 
by the thermally unstable, which moves outward from the inner to outer disk, perturbs the 
outer material which can produce the soft X-ray lag at a timescale of a few hundred seconds 
Li et al.l fl2007h . 
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Fig. 1. — The lightcurves (left), cross-correlations between the soft (2-3.3 keV) and hard (12- 
30 keV) X-rays (middle), the X-ray spectra (right) of the ObsID for hard and soft regions of 
the light curves are shown. 




Fig. 1. — Continued. 




Fig. 2. — CCD of all the analyzed observations is divided into four regions, where squares 
(dots, crosses) stand for the data with anti-correlation (ambiguous, positive correlation), 
respectively. 



Table 1: Parameters of ObsIDs where the anti-correlations are detected 



ObsID 


Date 
Year-Month-Day 


Location 


CC 


Delay(Error) 

(s) 


Hardness Ratio 


Energy of pivoting 
(keV) 


10068-06-01-00 


1996-09-01 


III/IV 


-0.34±0.02 


161±20 


4.69/0.60 


~4.5 


10068-06-02-00 


1996-09-04 


III/IV 


-0.37±0.01 


-350±10 


4.562/0.568 


~3.7 


20084-01-02-03 


1997-09-01 


IV 


-0.34±0.02 


-25.9±24.9 


4.636/0.592 


~6 


50029-34-02-00 


2000-08-14 


IV 


-0.29±0.03 


-16.8±21.0 


4.73/0.608 


~5 


70036-03-01-00 


2002-06-11 


IV 


-0.20±0.02 


317±8 


4.827/0.616 


~5.1 


91025-01-04-03 


2007-01-20 


IV 


-0.30±0.03 


88.7±11 


4.665/0.599 


~4.3 


91025-01-04-08 


2007-01-23 


III/IV 


-0.27±0.03 


63.4±33.4 


4.64/0.609 




91025-01-06-01 


2007-06-22 


III/IV 


-0.27±0.04 


-113±18 


4.571/0.589 


~4.2 


91025-01-06-02 


2007-06-23 


IV 


-0.27±0.03 


67.6±23.2 


4.61/0.607 


~5.1 


91025-01-06-08 


2007-06-27 


III/IV 


-0.25±0.02 


152±28 


4.584/0.601 


~3.9 


91025-01-06-09 


2007-06-28 


III/IV 


-0.30±0.02 


154±24 


4.495/0.593 


~4.7 


91025-01-09-05 


2008-01-13 


IV 


-0.18±0.02 


-192±49 


4.574/0.596 


~7 



Table 2: Percentage of positive and anti-correlation for each region of Fig. 1 





I 


II 


III 


IV 


positive-correlation 


19% 


90% 


94% 


28% 


anti-correlation 


0% 


0% 


2% 


30% 


ambiguous 


81% 


10% 


4% 


42% 



Table 3: The fitting parameters of the spectra of ObsID 91025-01-04-03 



Parameters 


fcT in (kcV) 


^disk 


kT bb 


N bb 


Flux diskbb / FluX tota i 


Flux bb /Flux total 


X 2 (d.o.f) 


a 
b 
c 
d 


1.95±0.09 
1.86±0.07 
1.73±0.05 
1.63±0.05 


18.6±2.8 
21.4±2.6 
25.7±2.7 
39.9±3.6 


2.72±0.05 
2.68±0.05 
2.65±0.04 
2.59±0.03 


0.051±0.003 
0.040±0.002 
0.031±0.002 
0.031±0.002 


4.07/8.23 (49.5%) 
3.91/7.22 (54.2%) 
3.32/5.91 (56.2%) 
2.75/5.51 (49.9%) 


4.12/8.23 (50.1%) 
3.26/7.22 (45.2%) 
2.55/5.91 (43.1%) 
2.66/5.51 (48.3%) 


1.20(45) 
0.84(45) 
1.20(45) 
1.31(45) 



The flux in unit of 10 9 ergs cm 2 s 1 is calculated in the energy band 2-30 keV. Errors are quoted at a 
90% confidence level. The letters a, b, c and d stand for the segments in the light curve. 
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Fig. 3.— Top: The light curves of soft and hard X-rays of ObsID 91025-01-04-03. Mid- 
dle and Bottom: The unfolded spectrum for each segment is plotted along with its model 
components. 



